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Abstract — A fundamental problem in wireless communica- 
tion is the time-frequency shift (TFS) problem: Find the time- 
frequency shift of a signal in a noisy environment. The shift is 
the result of time asynchronization of a sender with a receiver, 
and of non-zero speed of a sender with respect to a receiver. 
A classical solution of a discrete analog of the TFS problem 
is called the matched filter algorithm. It uses a pseudo-random 
waveform S(t) of the length p, and its arithemtic complexity 
is 0{p^ • log(p)), using fast Fourier transform. In these notes 
we introduce a novel approach of designing new waveforms that 
allow faster matched filter algorithm. We use techniques from 
group representation theory to design waveforms 5'(t), which 
enable us to introduce two fast matched filter (FMF) algorithms, 
called the flag algorithm, and the cross algorithm. These methods 
solve the TFS problem in 0(p • log(p)) operations. We discuss 
applications of the algorithms to mobile communication, GPS, 
and radar. 

I. Introduction 

Denote by H = C{¥p) the vector space of complex valued 
functions on the finite field ¥p = {0, 1, — 1}, where 
addition and multiplication is done modulo the odd prime 
number p. The vector space H is equipped with the standard 

inner product (/i,/2) = E fi{t)f2{t), for /i,/2 G and 

te¥p 

will be referred to as the Hilbert space of digital signals. 
Let us start with a motivational problem. 

A. Mobile communication problem 

We consider the following mathematical model of mobile 
communication ifTOl . There exists a collection of users j = 
l,...,r, each holding a bit 6 j G {±1}, and a private signal 
Sj eH. User j transmits its message bj • Sj to a base station 
(antenna), and the base station receives the superposition sum 

r 

R{t) = ^bj-e^^^-' -Sjit^Tj) +W(t), te¥p, (I.l) 

where W G H denotes a random white noise of mean zero, 
Tj encodes the time asynchronization of user j with the base 



station, ujj encodes the radial velocity of user j with respect 
to the base station, and i = 

The base station "knows" the signals ^^'s and R. The 
objective is: 

Problem I.l (Mobile communication problem): Extract the 
bits bj, j = l,...,r. 

A resolution of Problem |LT] will be deduced (see Section 
I-F) from our solution to the following problem. 



B. The time-frequency shift (TFS) problem 

We have r signals Sj G H, j = 1, ...,r, called the sender 
waveforms. Additionally, we are given the receiver waveform 
R eH, which satisfies 



R{t) 



e V ^ 



Sj{t^rj) +>V(t), te¥p, (1.2) 



where W G H denotes a random white noise of mean zero, and 
(rj^ujj) G ¥pX¥p, j = 1, r. We will call the pairs {rj^ujj) 
the time -frequency shifts, and the vector space V = ¥p x ¥p 
the time -frequency plane. 

The precise formulation of the time -frequency shift problem 
is the following: 

Problem 1.2 (TFS problem): Given the waveforms Sj, j = 
1, r, and R, extract the time-frequency shifts (rj^uoj) G V, 
j = l,...,r. 

C. The matched filter (MF) algorithm 

A classical solution O, El, Q, QOl, lEl, C3 to 
Problem [Oj is the matched filter algorithm. For a fixed k G 
{1, r}, we define the following matched filter (MF) matrix 
of the sender Sk, and the receiver R: 

M[Sk,R]{T,uj) = (e'l^'^-' • Suiter), R(t)), (t.uj) G V. 

(1.3) 



A direct verification shows that for 
1, r, we have 



e p 



M[Sk,R]{r,u) 



Cj^-M[Sk,Sk]{r-rk,uj-uJk) (1.4) 



■EC. 



M[Sk^Sj]{T 



+0(^1, 

where NSR = ^^r^ is the inverse of the signal-to-noise ratio 
between the waveform Sk and W. For simplicity, we assume 
that the NSR is not too large, and, for the rest of the paper, 
we will omit the last term in ( |I.4| ). 

In order to extract the time-frequency shift (t/cCJ/c), using 
the matched filter, it is "standard" (see 131. 141. iBIl. 171. ifTOl. 
ifTTl . (21) to use almost-orthogonal pseudo-random signals 
Sj eV, of norm one. Namely, all the summands in right-hand 



side of (1.4) are of size O(^), with the exception that for 

(r,cj) = (r/e,cjfe) we have M[Sk,Sk]{r -Tk.oJ- ojk) = 1. 
Hence, 



\M[Sk,R]{T,uj)\ 



/ 1 + ^r,p 



if (r,cj) = (r/e,cjfc); 
if (r,a;) ^ {Tk.uJk), 



where e. 



(1.5) 



the fastest method which exists in the literature IS). Note 
that computing one entry in M[SkjR] costs already 0{p) 
operations. This leads to the following fast matched filter 
(FMF) problem: 

Problem 1.4 (FMF problem): Design waveforms Sj G H, 
j = 1, ...,r, to solve TFS problem in almost linear time for 
shift. 

D. The flag method 

We introduce the flag method to propose a solution to FMF 
problem. We will show how to associate with the p + 1 lines, 
through (0, 0) in the time-frequency plane, Lj^ j = 1, 
a system of almost orthogonal waveforms Slj G that we 
will call flags. The system satisfies 

{2 + £^,p, if (r,a;) = (r/e,a;/e); 
l^Sr^p, if (r,cj) G L'^ \ (r/e,cj/e); 
Er^p, if (r,cj) G V \ L'^, 

(1.6) 

where Sr^p = O(^), is the receiver waveform (lO), defined 
with respect to any r flags containing Sl^i and is the 
shifted line + (rkjUJk)- 




Fig. 1. |Al[»S'i, i?] I with pseudo-random 81^82, and (ri,u:i) = (50,50) 

Identity ( [13] ) suggests the following "entry-by-entry" 
algorithmic solution to TFS problem: Compute the 
matrix M[Sk^R]^ and choose {rk^ojk) for which 
\M[Sk^R]{Tk^(jOk)\ ^ 1. However, this solution of TFS 
problem is very expensive in terms of arithmetic complexity, 
i.e., the number of arithmetic (multiplication, and addition) 
is 0{r • p^). One can do better using a "line-by-line" 
computation. This is due to the next observation. 

Remark 1.3 (FFT): The restriction of the matrix M[Sk, R] 
to any line (not necessarily through the origin) in the time- 
frequency plane is a convolution that can be computed, 
using the fast Fourier transform algorithm (FFT), in 0{p ■ 
log(p)) arithmetic operations. 

As a consequence of Remark Oj one can solve TFS 
problem in 0{r • p^ ■ \og{p)) arithmetic operations. To the 
best of our knowledge, the "line-by-line" computation is also 




Fig. 2. \M[Sl^,R]\ with two flags Sl^.Sl^, and (ri,a;i) = (50,50) 



Identity ( [L6| ) suggests the "flag" algorithmic solution to 
FMF problem in the case that the number of waveforms 
r <C ^/p^ and p is sufficiently large. In the following we 
assume that R and Sl^ are as in (O). 



Algorithm 1.5 (Flag algorithm): • Choose a line dif- 
ferent from Lfc. 

• Compute M[SLk-,R] on L^. Find (r,cj) such that 
\M[SLk-,R\{T.,uj)\ ^ 1, i.e., (r,cc;) on the shifted line 

• Compute M[SLkiR] on + {rk^ojk) and find (r,a;) 
such that \M[Sl,,,R]{t,uj)\ ^ 2. 

The arithmetic complexity of the flag algorithm is 0(r • 



plog(p)), using the FFT (Remark 1.3). 



E. The cross method 

Another solution to the TFS problem, and subsequently to 
the mobile communication problem, is the cross method. The 
idea is similar to the flag method, i.e., first to find a line on 




Fig. 3. Diagram of the flag algorithm 



which the time-frequency shift is located, and then to search 
on the Hne to find the time-frequency shift. We will show how 
to associate with the distinct pairs of lines L, M C ^ a 
system of almost-orthogonal waveforms S^^^^ that we will 
call crosses. The system satisfies 



2+£^,p, if (r,u;) = (r^,^,u;^,^); 

if (r,u;)G(L'UMO\(r^,^,^^,^, 
£r;p, if (r,u;)Gy\(L'UM'), 

where Sr^p = O(^), is the receiver waveform (Ol), defined 
with respect to any r different crosses containing^^^^ , and 




Fig. 4. \M[S^^,^^,R]\ with crosess S^^,^^,S^^ 
(50,50) 

The arithmetic complexity of the cross method is 0(r • 
plog(p)), using the FFT (Remark 1.3). 



E Solution to the mobile communication problem 

Looking back to Problem |lT] we see that the flag and 
cross algorithms suggest a fast 0{r ■ p ■ log(p)) solution to 
extract ALL the bits bk. Indeed, identity ( [L4| ) implies that 
bf. ^ [5'/c, i^](r/e, cj/e)/2, k = l,...,r, where R is the 
waveform (10), with Sj = Slj, j = 1,...,^, for the flag 
method, or^j = ,m ' J = l^...^r, for the cross method, 
and r <C ^/p. 



of a pseudorandom signal and a structural signal. The first 
has the MF matrix which is almost delta function at the 
origin, and the MF matrix of the second is supported on a 
line. The designs of these waveforms are done using group 
representation theory. The pseudorandom signals are designed 
141 . O, 1 12] using the Weil representation, and will be called 
Weil (peak) signal^ The structural signals are designeded 
161 . fTl using the Heisenberg representation, and will be 
called Heisenberg (lines) signals. We will call the collection 
of all flag waveforms, the Heisenberg-Weil flag system. In 
this section we briefly recall constructions, and properties of 
these waveforms. A more comprehensive treatment, including 
proofs, will appear in |2|. 

A. The Heisenberg (lines) system 

Consider the following collection of unitary operators, 
called Heisenberg operators, that act on the Hilbert space of 
digital signals: 



7r(r,a;) M,^ o L^, 



(II.l) 



where L,[/](t) 



f{t + r) is the time-shift operator, 
f{t) is the frequency- shift operator, for 
every f e H, t e ¥p, and o denotes composition of operators. 



-UJ-t 



The operators (II.l) do not commute in general, but 
rather obey the Heisenberg commutation relations 7r(r,co') o 
7r(r',cc;') = e^''^^ ^ • 7r(r', cj') o 7r(r, a;). The expression 
Tuj' — ujt' vanishes if (r,cc;), {r' ^uo') belong to the same 
line. Hence, for a given line L C F = Fp x Fp we have 
a commutative collection of unitary operators 



(II.2) 



We use the theorem from linear algebra about simultaneous 
diagonalization of commuting unitary operators, and obtain 
O, Q a natural orthonormal basis Bl C H consisting of 



common eigenfunctions for all the operators (II.2). The system 
of all such bases Sl, where L runs over all lines through the 
origin in will be called the Heisenberg (lines) system. We 
will need the following result 161, (TJ: 

Theorem U.l: The Heisenberg system satisfies the proper- 
ties 

1) Line. For every line L (Z V , and every G Bl, we 
have 



1, if (t,uj) e L; 
0, if {r,uj) ^ L. 



2) Almost- orthogonality. For every two lines Li ^ L2 C 
V, and every Jl^ e Bl^, Jl^ G Bl^, we have 

1 



II. The Heisenberg-Weil flag system 



for every (r, a;) G V. 



The flag waveforms, that play the main role in the flag 
algorithm, are of a special form. Each of them is a sum 



^For the purpose of the Flag method, other pseudorandom signals may 
work. 




2) Almost- orthogonality. For every two tori Ti, T2 C G, 
and every ifrj._^ G Bt^ , (/^^^ ^ 7 with (/p^^ 7^ (/^^^ 5 we 
have 



for every (r, a;) G V. 



Fig. 5. for L = {(t,t); r € Fp} 



B. rfte Wfei7 (peaks) system 

Consider the following collection of matrices 

d I 



G = SL2{Wp) 



a, 6, c, d G Fp, and ad — be = 1 




Note that G is in a natural way a gr(9w/7 |1| with respect to 
the operation of matrix multiplication. It is called the special 
linear group of order two over Fp. Each element g ^ G acts 
on the time-frequency plane V via the change of coordinates 
V ^ g ■ V. For every g G G, let p{g) be a linear operator on 
1-L which is a solution of the following system of linear 
equations: 

T.g\ p{g)o7r{r,uj) = 7r{g-{r,uj))op{g)^ r, a; G F^, (II.3) 

where tt is defined by ( [II. Denote by Sol(E^) the space of 
all solutions to System pO] ). The following is a basic result 
CI: 

Theorem II. 2 (Stone-von Neumann-Schur-Weil): There ex- 
ist a unique collection of solutions {p{g) G Sol(i;^); g e G}^ 
which are unitary operators, and satisfy the homomorphism 
condition p{g • h) = p{g) o p{h). 

Denote by U{l-L) the collection of all unitary operators 



on the Hilbert space of digital signals H. Theorem |II.2 
establishes the map p : G U{l-L), which is called the 
Weil representation |[T3]| . The group G is not commutative, but 
contains a special class of maximal commutative subgroups 
called tor|3 PI, PI. Each torus T c G acts via the Weil 
representation operators 



p{g) :n^n, geT. 



(11.4) 



This is a conmiutative collection of diagonalizable operators, 
and it admits f4l, f5| a natural orthonormal basis Bt for H, 
consisting of common eigenf unctions. The system of all such 
bases St, where T runs over all tori in G, will be called the 
Weil (peaks) system. We will need the following result |4J, ll5J: 
Theorem II. 3: The Weil system satisfies the properties 
1) Peak. For every torus T C G, and every (pj. g Bt, we 
have 

1, if (r,a;) = (0,0); 
if (r,a;)/(0,0). 



|A^[(^^,(^^](r,cj)| 



Fig. 6. M[^t.^t]^otT = {(^^ ^^i);0/aGFp} 



C. The Heisenherg-Weil system 

We define the Heisenberg-Weil system of waveforms. This 
is the collection of signals in which are of the form 
Sl = fh^^T^ where and (/p^ are Heisenberg and Weil 
waveforms, respectively. The main technical result of this 
paper is: 

Theorem II. 4: The Heisenberg-Weil system satisfies the 
properties 

1) Flag. For every line L C torus T C G, and every 
flag Sl = fh^ with Jl G Sl, ^ ^t, we have 



\M[Sl^Sl\{t,uj)\ 



if (r,a;) = (0,0); 
if (r,cj) G L \ (0,0); 



2 + 
1 + 

Sp, if (r,cj) G F \ I/, 



where \ep\ < ^, and \ep\ < ^. 
2) Almost- orthogonality. For every two lines Li ^ L2 C 
tori Ti, T2 C G, and every two flags Sl^ = /l, + 

V^T,, with /l^. G Sl,-, ^t, ^ >St,., J = 1, 2, (/:)^^ 7^ 
we have 



\M[Sl,^Sl,]{t, 




ifTi/T2; 
ifTi=T2, 



for every (r, a;) G V. 
A proof of Theorem |II.4| will appear in O . 



^ There are order of tori in SL2{¥p). 



Remark II. 5: As a consequence of Theorem II.4 we obtain 
families of p + 1 almost-orthogonal flag waveforms which 
can be used for solving the TFS and mobile communication 
problems in almost linear time. 




Fig. 7. \M[Sl,Sl]\ for Heisenberg-Weil flag with L = {(r, r); r G Fp} 



III. The Heisenberg cross system 

We define the Heisenberg cross system of waveforms. This 
is the collection of signals in H, which are of the form S*^ ^ = 
/l + /m, where /l, /m, L ^ are Heisenberg waveforms 
defined in Section |II-A[ The following follows immediately 
from Theorem III. II 

Theorem III.l: The Heisenberg cross system satisfies the 
properties 

1) Cross. For every pair of distinct lines L, M C V, and 
every cross 6'^,^ Il^Im, with /l G Sl, /m ^ Sm, 
we have 



where < 



2+£p, if (r,a;) = (0,0); 

if (r,u;)G(LUM)\(0,0); 
£p, if (r,a;)Gy\(LUM), 



to the earth its signal Sj G H multiplied by its location hj. 
The client receives the signal 

r 

where ujj encodes the radial velocity of satellite j with respect 
to the client, Tj encodes the distance between satellite j and 
the clien|5 and W is a random white noise of mean zero. 

Problem IV. 1 (GPS problem): Find {bj^rj)^ j = 1, ...,r. 

By using Heisenberg-Weil or Heisenberg cross waveforms 
we find the pairs {bj^rj) in 0(r • plog(p)) arithmetical 
operations. 

B. Application to discrete radar 

The model of discrete radar works as follows |T|. A radar 
sends a waveform S e H which bounds back by r targets. 
The signal R e 1-L which is received as an echo has the forn]^ 

r 

where ujj encodes the radial velocity of target j with respect 
to the radar, tj encodes the distance between target j and the 
radar, and W is a random white noise of mean zero. 

Problem IV.2 (Discrete radar problem): Find (rj^ujj)^ j = 
l,...,r. 

By sending Heisenberg-Weil waveform Sl = fh ^ 



2) Almost- orthogonality. For every four distinct lines 

Li,Mi,L2,M2 C V, and every two crosses ,^ = 
/l, + /m, , j = 1, 2, we have 



< 



for every (r, cj) G V. 



Remark III. 2: As a consequence of Theorem III.l we obtain 
families of almost-orthogonal cross waveforms which 
can be used for solving the TFS and mobile communication 
problems in almost linear time. 

IV. Applications to GPS and radar 

In the introduction we described application of flag and 
cross methods to mobile communication. In this section we 
demonstrate applications to global positioning system (GPS), 
and discrete radar. 

A. Application to global positioning system (GPS) 

The model of GPS works as follows fS^. A client on the 
earth surface wants to know his geographical location. Satel- 
lites j = 1 , . . . , r send to earth their location. For simplicity, the 
location of satellite j is a bit 6j G {±1}. Satellite j transmits 



geff 



\M[SL,R]{r,co)\ 



2 + if (r,cj) G {{Tj.Uj)}; 
1 + er,p, if (r,cj) G L + {rj.cuj) \ {rj.ujj)', 
£rv^ otherwise. 



where e. 



This means that by using the flag algorithm we solve the 
radar problem in 0{r ■ p\og{p)) arithmetical operations. 

Remark IV. 3 (Important): Note that the cross method is not 
applicable for the discrete radar problem if the number of 
targets r > 1. 
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^From the tj we can find 1 8 1 the distance between satellite j and the client, 
given that r > 4 and the clocks of all satellites are synchronized. 

^In practice there are intensity coefficients < < 1 such that R{t) = 

aye^^ -Sit + Tj) +W(t). Assuming that aj 's are sufficiently 
large our methods are applicable verbatim. 

^For simplicity we assume that all the shifted lines L + (rj,cjj)'s are 
distinct. The general case is treated similarly. 




Fig. 8. \M[Sl,R]\ with L = {(r,0); r G Fp}, and shifts (50,50), 
(100,100), (150,150) 
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